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Abstract 

Let £ = {x^}"^ be a grid of n points in the d-cube I d :— [0, l] d , and $ = {Vj}?=i a family 
of n functions on l d . We define the linear sampling algorithm L„($,£, •) for an approximate 
recovery of a continuous function / on I d from the sampled values /(a: 1 ), f(x n ), by 

n 
3=1 

For the Besov class e of mixed smoothness a (defined as the unit ball of the Besov space 
MB" ( ), to study optimality of L n (<I>,£, •) in £ 9 (I d ) we use the quantity 

r„(B« e ) g := inf sup [|/- £«(#,£,/)||„ 

where the infimum is taken over all grids £ = {x-?}™ =1 and all families $ = {^}" = i in L g (I d ). 
We explicitly constructed linear sampling algorithms L„(<£>,£, •) on the grid £ = G d (m) := 
{(2~ fcl si, 2~ kd Sd) 6 I d : fci + ... + kd < m}, with $ a family of linear combinations of mixed 
B-splines which are mixed tensor products of either integer or half integer translated dilations of 
the centered B-spline of order r. The grid G d (m) is of the size 2 m m d_1 and sparse in comparing 
with the generating dyadic coordinate cube grid of the size 2 dm . For various < p, q, 9 < oo 
and 1/p < a < r, we proved upper bounds for the worst case error supj eS c« e ||/ — L n (<I>,£, f)\\ q 
which coincide with the asymptotic order of r n (B" e ) q in some cases. A key role in constructing 
these linear sampling algorithms, plays a quasi-interpolant representation of functions / G B" e 
by mixed B-spline series with the coefficient functionals which are explicitly constructed as linear 
combinations of an absolute constant number of values of functions. Moreover, we proved that 
the quasi- norm of the Besov space MB^ e is equivalent to a discrete quasi- norm in terms of the 
coefficient functionals. 

Keywords Linear sampling algorithm • Quasi-interpolant • Quasi-interpolant representation • 
Mixed B-splinc • Besov space of mixed smoothness. 

Mathematics Subject Classifications (2000) 41A15 • 41A05 • 41A25 • 41A58 • 41A63. 



1 Introduction 



The aim of the present paper is to investigate linear sampling algorithms for recovery of functions 
on the unit d-cube I d := [0, l] d , having a mixed smoothness. Let £ = {x- J }™ =1 be a grid of n points 
in I d , and <3? = {<£;}? = i a family of n functions on I d . Then for a continuous function / on I d , we 
can define the linear sampling algorithm L n = L n ($>,£, •) for approximate recovering / from the 
sampled values f(x 1 ),...,f(x n ), by 

n 

L n (f) = L n ($,t,f):= Y,f( xJ )Vj- (1.1) 

i=i 

Let L q := L q (I d ), < q < oo, denote the quasi-normed space of functions on I d with the gth 
integral quasi- norm || • \\ q for < q < oo, and the ess sup- norm || • for (7 = 00. The recovery 
error will be measured by ||/ — L n (<f>,£, f)\\ q . 

If W is a class of continuous functions, sup^g^a^ \\f — L n ($,£, f)\\ q is the worst case error of the 

recovery of functions / from W by the linear sampling algorithm L n (3>,£, •). To study optimality 
of linear sampling algorithms of the form (jl.ip for recovering / £ W from n their values, we will 
use the quantity 

r n (W) q := inf sup \\f - L n {*,Z,f)\\ q , (1.2) 

where the infimum is taken over all grids £ = {x : '}™ =1 and all families = {^j}" =1 in L q . 

A challenging problem in linear sampling recovery of functions from a class W with a given 
mixed smoothness, is to construct a sampling algorithm L ra (<3?,£, •) with an appropriate sampling 
grid £ = {x : '}™ =1 and family <3? = {^j}™ =1 which would be asymptotically optimal in terms of the 
quantity r n (W) q . 

For periodic functions Smolyak [23J first constructed a specific linear sampling algorithm based 
on the de la Vallee Poussin kernel and the following dyadic grid in I d 

G d (m):= {(2- fcl si,...,2- fed s d ) el d : k E A(m)} = {2~ k s : k G A(m), s G I d (k)}. 

Here and in what follows, we use the notations: xy := (xiy\, ...,Xdl/d)', 2 X '■= (2 X1 , ...,2 Xd ); \x\i : = 
Y$=i \xi\ for x,y £ R d ; A(m) := {ic£Z|: < m}; I d {k) := {s £ Z d + : < Si < 2 k \ i £ N[d]}; 
N[d] denotes the set of all natural numbers from 1 to d; X{ denotes the ith coordinate of x £ M. d , 
i.e., x := (xi, Xrf). Temlyakov [23], [26], [27] and Dinh Dung [9j-[TT] developed Smolyak's 
construction for study the asymptotic order of r n (W) q for periodic Sobolev classes and Holder 
classes as well their intersection. In particular, the first asymptotic order 

r n {H*) q x (n- 1 log d - 1 n) a - 1 / p+1 / q (log d - 1 n)V q , K P <q<2, a>l/p, 

was obtained in [9]-|10j. For non-periodic functions of mixed smoothness 1/p < a < 2, this 
problem has been recently studied by Sickel and Ullrich [22j, using the mixed tensor product of 
piecewise linear B-splines (of order 2) and the grid G d (m). It is interesting to notice that the linear 
sampling algorithms considered by above mentioned authors are interpolating at the grid G d (m). 
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Naturally, the quantity r n (W) q of optimal linear sampling recovery is related to the problem 
of optimal linear approximation in terms of the linear n-width X n (W) q introduced by Tikhomirov 
[28]: 

\ n {W) q := inf sup \\f - A n (f)\\ q , 

where the infimum is taken over all linear operators A n of rank n'mL q . The linear n-width \ n (W) q 
was studied in [14], [20], [21], ect. for various classes W of functions with mixed smoothness. The 
inequality r n > A n is quite useful in investigation of the (asymptotic) optimality of a given linear 
sampling algorithm. It also allows to establish a lower bound of r n via a known lower bound of 

In the present paper, we continue to research this problem. We will take functions to be 
recovered from the Besov class 5" g of functions on I d , which is defined as the unit ball of the 
Besov space MB^g having mixed smoothness a. For functions in B^g, we will construct linear 
sampling algorithms L n (<I>,£,-) on the grid £ = G d (m) with <I> a family of linear combinations 
of mixed B-splines which are mixed tensor products of either integer or half integer translated 
dilations of the centered B-spline of order r > a. We will be concerned with the worst case error 
of the recovery of B^ e in the space L q by these linear sampling algorithms and their asymptotic 
optimality in terms of the quantity r n (Bp g ) q for various < p, q, 6 < oo and 1/p < a < r. A key 
role in constructing these linear sampling algorithms, plays a quasi-interpolant representation of 
functions / £ MB^g by mixed B-spline series which will be explicitly constructed. Let us give a 
sketch of the main results of the present paper. 

We first describe representations by mixed B-spline series constructed on the basic of quasi- 
interpolants. For a given natural number r, let M be the centered B-spline of order r with support 
[— r/2, r/2] and knots at the points —r/2, — r/2+1, r/2 — 1, r/2,. We define the integer translated 
dilation M kjS of M by 

M h ,s(x) := M(2 k x - s), k £ Z + , sGZ, 
and the mixed d-variable B-spline Mf. s by 

d 

M ktS {x) := Y[M ki , Si { Xi ), kez d + , se Z d , (1.3) 

i=l 

where Z + is the set of all non-negative integers, := {s £ Z d : Si > 0, i £ N{d}}. Further, we 
define the half integer translated dilation of M by 

Ml s (x) := M(2 k x - a/2), k £ Z + , s £ Z, 
and the mixed d- variable B-spline by 

d 

M fc * s (x) := \{M^ Si {x t ), keZ d + , S £ Z d . 
i=i 

In what follows, the B-spline M will be fixed. We will denote M£ := Mfc iS if the order r of M is 
even, and M r k s := M£ g if the order r of M is odd. 
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Let < p, 6 < oo, and 1/p < a < min(r, r — 1 + 1/p). Then we prove the following mixed 
B-spline quasi-interpolant representation of functions / G MB^ e . Namely, a function / in the 
Besov space MB^ e can be represented by the mixed B-spline series 

/ = E E m 

converging in the quasi-norm of MB^g, where J d (k) is the set of s for which M£ s do not vanish 
identically on I d , and the coefficient functionals el (/) explicitly constructed as linear combinations 
of an absolute constant number of values of / which does not depend on neither k, s nor /. More- 
over, we prove that the quasi-norm of MB^ g is equivalent to some discrete quasi-norm in terms 
of the coefficient functionals c£ .(/). B-spline quasi-interpolant representations of functions from 
the isotropic Besov sapces has been constructed in [12] . |13j . Different B-spline quasi-interpolant 
representations were considered in [7]. Both these representations were constructed on the ba- 
sic of B-spline quasi-intepolants. The reader can see the books [2], [5] for survey and details on 
quasi-interpolants. 

Let us construct linear sampling algorithms L n (<]?,£, •) on the grid £ = G d (m) on the basic of 
the representation (II. 4|) . For m € Z+, let the linear operator R m be defined for functions / on I d 
by 

«"»(/) = E E c lM) M *,s- (i.5) 

fceA(m) sGJ^(fe) 

If fh is the largest of m such that 

2 m m d-i x |G*(th)| < n 

for a given n, where \A\ denotes the cardinality of A, then the operator is a linear sampling 
algorithm of the form (jl.ip on the grid G d {m). More precisely, 

= L n (<D,£,/) = ^ f(2~ k s)^ s , 

where are explicitly constructed as linear combinations of an absolute constant of B-splines 
M£ ., which does not depend on neither k, s nor /. It is worth to emphasize that the grid G d {m) 
is of the size 2 m m d ~ 1 and sparse in comparing with the generating dyadic coordinate cube grid 
of the size 2 dm . We give now a brief of our results concerning with the worst case error of the 
recovery of functions / from B® e by the linear sampling algorithms Rm(f) and their asymptotic 
optimality. 

We use the notations: x + := max(0,x) for x G R; A n (f) < B n (f) if A n (f) < CB n (f) with C 
an absolute constant not depending on n and/or / € W, and A n (f) x B n (f) if ^4 n (/) B n (f) 
and B n (f) <^ A n (f). Let us introduce the abbreviations: 

E(m):= sup ||/ - i? m (/) ||„ r n := r n (B" 9 ),. 

Let < p,q,8 < oo and 1/p < a < r. Then we have the following upper bound of r n and 
£?(m). 
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(i) For p > q, 

r n < E(m) < { 

(ii) For p < q, 

r n < E(m) < 



(n- 1 log d ~ 1 n) a , 
(n- 1 log d - l n) a (log d - 1 
(7^- 1 log <^ - 1 r^) a (log , 



n 



,1/5-1/ 



d-i n )i-i/e 



< min(g, 1), 

> min((7, 1), q < 1, 

> min((7, 1), g > 1. 



(1.6) 



' (n- 1 log^ 1 n) a - 1 IP+ 1 li(\og d - x n)( 1 /<?-i/e)+ ; 



(n 1 loj 



.d-l 



n 



.*-i n )(i-i/0) H 



g < oo, 
g = oo. 



(1.7) 



From the embedding of MB^ g into the isotropic Besov space of smootthness da and known 
asymptotic order of the quantity (jl.2p of its unit ball in L q (see [5] , [16] , [17] , [18] , [27] ) it follows 
that for < p, q < oo, < 9 < oo and a > 1/p, there always holds the lower bound r n 
n -a+(i/p-i/ij) + _ However, this estimation is too rough and does not lead to the asymptotic order. 
By use of the inequality Xn(Bp d ) q > r n and known results on \ n {Bpo)q [15] . [20] . from fjl .6|) and 
(jl.7p we obtain the asymptotic order of r n for some cases. More precisely, we have the following 
asymptotic orders of r n and E(fh) which show the asymptotic optimality of the linear sampling 
algorithms Rm- 



(i) For p > q and 9 < 1, 

E(fh) x r n x (rT L log a_1 n 

(ii) For 1 < p < q < oo, 



-ii_d-l„v* \2<q<P<°o, 



1 < p = q < oo. 



(1- 



£(m) x r n x (7 l - 1 log d - 1 r i ) a - 1 /f+V 9 (i og rf-i n )(i/ g -Ve) +) J 2 ^< 2 < < g, ( 

[9 < 2. 

The present paper is organized as follows. In Section [2] we give a necessary background of 
Besov spaces of mixed smoothness, B-spline quasi-interpolants, and prove a theorem on the mixed 
B-spline quasi-iterpolant representation (jl.4p and a relevant discrete equivalent quasi-norm for the 
Besov space of mixed smoothness MB^ e . In Section [3] we prove the upper bounds (ll.6j) - (ll.7j) and 
the asymptotic orders f)1.8|) — f)1.9[) . In Section [JJ we consider interpolant representations by series 
of the mixed tensor product of piecewise constant or piecewise linear B-splines. In Section [5] we 
present some auxiliary results. 



2 B-spline quasi-interpolant representations 

Let us introduce Besov spaces of functions with mixed smoothness and give necessary knowledge 
of them. For univariate functions the Zth difference operator is defined by 

Km--= £(-i)'- j Y'.W+j>o- 

j=0 VJ/ 
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If e is any subset of N[d], for multivariate functions the mixed (7,e)th difference operator A^ e is 
defined by 

where the univariate operator is applied to the univariate function / by considering / as a 
function of variable Xi with the other variables held fixed. For a domain in W 1 , denote by L p {Vt) 
the quasi-normed space of functions on fi with the pth. integral quasi-norm || • || Pj q for < p < oo, 
and the ess sup-norm || • ||oo,f2 for p = oo. Let 

wf(/,t) p := sup ||A^(/)|| PiId(M , tel d , 

\hi\<ti,iee 

be the mixed (/, e)th modulus of smoothness of /, where I d (h, e) := {x € I d : Xi, X{ + lhi 6 I, i € e} 
(in particular, uf(f,t) p = ||/|| p ). We will need the following modified (Z, e)th mixed modulus of 
smoothness 

i/p 



Viet Ju(t,e)Jl"(h,e) , 



where U(t,e) := {x £l d : \xi\ <t, i € e}. There hold the following inequalities 

C lW f(f,t) p <^ e (/,t) P < C 2 w?(f,t) p (2.1) 

with constants Ci, C2 which depend on Z,p, d only. A proof of these inequalities for the univariate 
modulus of smoothness is given in [19] . They can be proven in a similar way for the multivariate 
(I, e)th mixed modulus of smoothness. 

If < p,9 < 00, a > and I > a, we introduce the quasi-semi-norm |/|n a > e for functions 
f S L p by 



(in particular, |/j 



MB a, S ~ IIJ IIP-*- 



For < p, 9 < 00 and < a < I, the Besov space MB® e is defined as the set of functions 
/ £ L p for which the Besov quasi-norm ||/|[m,b q s is finite. The Besov quasi-norm is defined by 

B{f) = H/IImb", := Yl \f\MB«>°- 

eCN[d\ 

We will study the linear sampling recovery of functions from the Besov class 

Kfi ■= if G MB P,e ■ B (f) < !>> 

with the restriction on the smoothness a > 1/p, which provides the compact embedding of MB^ e 
into C(I d ), the space of continuous functions on I d with max-norm. We will also study this 
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problem for B£ e with the restrictions a = 1/p and p < min(l, 9) which is a sufficient condition for 
the continuous embedding of MB^ d into C(I d ). In both these cases, B^ g can be considered as a 
subset in C{l d ). 

For any e C N[d], put Z£(e) := {s E Z^„ : s; = 0, i £ e} (in particular, Z^(0) = {0} 
and Z^(iV[d]) = Z 1 ^). If {gk}k£Z'l (e) ^ s a sequence whose component functions ^ are in L p , for 

< p, 9 < oo and /3 > we define the bg' e (L p ) "quasi-norms" 

Hfe}ll r(ip) - ( E ( 2/3|fc|l lb fc || P ) 
v fcez;[(e) 

with the usual change to a supremum when 9 = oo. When {gk}kez d ( e ) ^ s a positive sequence, we 
replace \\gk\\p by \gk\ and denote the corresponding quasi- norm by ||{gjfe}||^,e. 

For the Besov space MB^ g , from the definition and properties of the mixed (Z,e)th modulus 
of smoothness it is easy to verify that there is the following quasi-norm equivalence 

B(f) x Bt(f) := £ ||K(/,2-*) p }|| 6r . 

eCN[d] 

Let A = {A(s)} Jg p( /t ) be a finite even sequence, i.e., A(— j) = X(j), where P((J>) 
\j\ < /it} and /i > r/2 — 1. We define the linear operator Q f° r functions / on R by 

:= ^A(/, S )M(x-s), 

where 

A(/, S ):= 2 

The operator Q is bounded in C(R) and 

||Q(/)||c(R) < l|A|lll/lb(R) 

for each / E C(R), where 

l|A||= E 

Moreover, Q is local in the following sense. There is a positive number 5 > such that for any 
/ G C(R) and x € R, Q(f,x) depends only on the value f(y) at an absolute constant number of 
points y with \y — x\ < 5. We will require Q to reproduce the space V r -i of polynomials of order 
at most r — 1 , that is, Q(p) = p, p € V r -i- An operator Q of the form (|2.2p - (j2.3p reproducing 
V r -ii is called a quasi-interpolant in C(R). 

There are many ways to construct quasi-interpolants. A method of construction via Neumann 
series was suggested by Chui and Diamond [3] (see also [2J p. 100-109]). A necessary and sufficient 



:= {j € Z : 
(2.2) 

(2.3) 
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condition of reproducing V T -\ for operators Q of the form (|2.2p - (|2.3p with even r and > r/2, was 
established in pQ. De Bore and Fix [4] introduced another quasi-interpolant based on the values 
of derivatives. 

Let us give some examples of quasi-interpolants. The simplest example is a piecewise constant 
quasi-interpolant which is defined for r = 1 by 

Q(f,x) := J2f(s)M(x-s), 

sez 

where M is the symmetric piecewise constant B-spline with support [—1/2, 1/2] and knots at the 
half integer points —1/2, 1/2. A piecewise linear quasi-interpolant is defined for r = 2 by 

Q(f,x):= £/( s )M(x-a), (2.4) 

where M is the symmetric piecewise linear B-spline with support [—1, 1] and knots at the integer 
points —1,0,1. This quasi-interpolant is also called nodal and directly related to the classical 
Faber-Schauder basic. We will revisit it in Section UJ A quadric quasi-interpolant is defined for 
r = 3 by 

Q(f, x) := \{~fi s ~ !) + " f( s + !) " 

where M is the symmetric quadric B-spline with support [—3/2, 3/2] and knots at the half integer 
points —3/2, —1/2, 1/2,3/2. Another example is the cubic quasi-interpolant defined for r = 4 by 

where M is the symmetric cubic B-spline with support [—2, 2] and knots at the integer points 
-2,-1,0,1,2. 

If Q is a quasi-interpolant of the form fj2.2[) — fj2.3|) . for h > and a function / on R, we define 
the operator Q(-;h) by 

Q{f;h) = a h oQ oa 1/h (f), 
where Uh{f,x) = f(x/h). By definition it is easy to see that 

Q(f,x;h)= Y,Mf,k;h)M(h- l x-k), 
k 

where 

A(f,k;h):= Y X (j)f(Hk-j)). 

The operator Q(-;h) has the same properties as Q: it is a local bounded linear operator in 
C(R) and reproduces the polynomials from V r -i- Moreover, it gives a good approximation for 
smooth functions [5, p. 63-65]. We will also call it a quasi-interpolant for C(R). However, the 
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quasi-interpolant Q(-;h) is not defined for a function / on I, and therefore, not appropriate for an 
approximate sampling recovery of / from its sampled values at points in I. An approach to con- 
struct a quasi-interpolant for functions on I is to extend it by interpolation Lagrange polynomials. 
This approach has been proposed in [12] for the univariate case. Let us recall it. 

For a non-negative integer k, we put Xj = j2 ,j € Z. If / is a function on I, let 

r-1 2 sk A s f(xn) S_1 
U k (f,x):= /(so) + £ 2^1111 II^-^)' 

r-i 2 sk A s k f(x 2 k r i) s ~ l 
V k (f,x):= f{x 2k _ r+1 ) + 2 " fc , 2 ~ r+1 n( x_x 2*-r+i+i) 

8=1 S ' j=0 

be the (r — l)th Lag range polynomials interpolating f at the r left end. points xo ; xi ; and r 
right end points x 2 fe -r+i> x 2 fc -r+3> ...,x 2 fc, of the interval I, respectively. The function f k is defined 
as an extension of / on R by the formula 



fk(x) :-- 



U k {f,x), x<0, 
fix), 0<x<l, 

y k (f,x), x>i. 



Obviously, if / is continuous on I, then f k is a continuous function on R. Let Q be a quasi- 
interpolant of the form (pT2"j) - (pT3|) in C(R). Put Z + := {fc € Z : k > -1}. If A; € Z+, we introduce 
the operator Q k by 

Qk(f,x)= Q(f k ,x;2- k ), and Q-i(f,x) := 0, x E I, 
for a function / on I. We have for € Z + , 

Qk(f,x) = Yl a M (/)M M (x), VxGl, (2.5) 

where </(£;) := {s € Z : — r/2 < s < 2 fc + r/2} is the set of s for which -Mfc jS do not vanish 
identically on I, and the coefficient functional a/^s is defined by 

o M (/):= A(A, S ;2" fc )= £ A(j)/ fc (2- fc ( 5 - ;)). 

lil<M 

Put Z+ := {fc £ Z{ : fej > — 1, ieN[d]}. For fc G Z^, let the mixed operator Q fc be defined by 

<2 



Qk — YlQkt, (2.6) 



1=1 



where the univariate operator Q ki is applied to the univariate function / by considering / as a 
function of variable Xi with the other variables held fixed. 
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We have 

Qk(M = Yl a kAf) M kA x )> Vx € i d , 

s£j d (k) 

where M k , s is the mixed B-spline defined in (|1.3p . J d (k) := {s £ Z d : — r/2 < Sj < 2 ki +r/2, i G 
iV[<i]} is the set of s for which M/% jS do not vanish identically on I d , 

ak,s(f) '■= a fci, S i((afe 2 , S2 (-afe d , Sd (/))), (2.7) 

and the univariate coefficient functional afc i)Si is applied to the univariate function / by considering 
/ as a function of variable Xi with the other variables held fixed. 

The operator Qk is a local bounded linear mapping in C(I d ) and reproducing V d _ 1 the space of 
polynomials of order at most r — 1 in each variable X{. More precisely, there is a positive number 
5 > such that for any / G C(I d ) and x G I d , Q(f,x) depends only on the value f(y) at an 
absolute constant number of points y with |yj — X{\ < 52~ ki ,i G iV[d]; 

||Qfc(/)||c(^)<C||A|| d ||/|b(i d ) (2.8) 

for each / G C^F*) with a constant C not depending on A;; and, 

Qkip*) = P,P^V d _i, (2.9) 

where p* is the restriction of p on The multivariate Qk is called a mixed quasi-interpolant in 
C(I d ). 

From (|2.8p and (j2.9j) we can see that 

||/-Q fc (/)|| C(Id ) ^0, fc^oo. (2.10) 

(Here and in what follows, k —> oo means that &j — > oo for z G -/V"[c?]). 

If /c G Z+, we define Tk := I — Qk for the univariate operator Q^, where I is the identity 
operator, if k G Z+, we define the mixed operator Tk in the manner of the definition (|2.6p by 

i=l 

For any e C N[d], put Z+(e) := {s G Z^ : Sj > —1, i G e, Sj = —1, i ^ e} (in particular, 
Z|(0) = {(-1, -1, -1)} and Z^(iV[cZ]) = Z£). We have Z^it) n Z%(v) = if u ± v, and the 
following decomposition of Z 6 ^: 

Z^ = (J Z*(e). 

If r is a number such that < r < min(p, 1), then for any sequence of functions {gk} there is 
the inequality 

||J>J T < El^llp- ( 2 - n ) 
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Lemma 2.1 Let < p < oo and r < min(p, 1). XTien /or any f € C(I rf ) and k G Z+(e), i/iere 
/io/ds i/ie inequality 



\T k (f)\\ p < C 



( \ 

^seZ^(e), s>k i 



1/r 



(2.12) 



TOt/i some constant C depending at most on r,^,p,d and ||A||, whenever the sum in the right-hand 
side is finite. 

Proof. Notice that Zi(0) = {(—1, —1, —1)} and consequently, the inequality ()2. 12|) is trivial for 
e = 0: < Cujf(f,l)p = C||/|| p . Consider the case where e ^ 0. For simplicity we prove 

the lemma for d = 2 and e = {1,2}, i.e., Z+(e) = Z?_. This lemma has proven in |12| . |13j for 
univariate functions (d = 1) and even r. It can be proven for univariate functions and odd r in a 
completely similar way. Therefore, by f)2. 1 1) there holds the inequality 



\T kl {f))\\ P « ( E I 2{Si ~ ki),p \ 2 ~ Si [ I K(/>*)i p ^ dh i 



, z = 1,2, 



, (2 ' 13) 

where the norm 117^(7)) || p is applied to the univariate function / by considering / as a function 
of variable Xj with the other variable held fixed. 

If 1 < p < oo, we have by (|2.13p applied for i = 1, 

VpI p 
il ? dx 2 



\T kl T k2 (f))\\ P « I [I E 2^"^ [2-* / / \A l hl ((T'f),xW dx x dh 



2( s i _A; i)/p ( 2 _Sl / 



l hl ((Tij),x)\P dx, dh dx 2 ) 

)./:(/,,, - y 

(T fc2 (A^/),x)| p dx 2 ) dh^\ . 



)-/n(/n) Ul 



Hence, applying (I2.13P with i = 2 gives 

f \A l h (f,x)\v dxdh) 

U{2- a )Jl 2 (h) I 



\T kl T k2 (f))\\ p « ^2l s - fc l^(2-l^ / 

< ^2l s - fe l 1 )/ p 'u; r (/,2-' £ ) p < E 2l '~ fc|l)/Pw »-(/> 2 ~*) 



v 

s>k s>k 



Thus, the lemma has proven when 1 < p < oo. The cases < p < 1 and p = oo can be proven in 
a similar way. □ 
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Let J d (k) := J d (k) if r is even, and J r d (/fc) := {s e Z d : -r < s, < 2 fe ' +1 +r,iG JV[d]} if r is 
odd. Notice that J d (k) is the set of s for which MI do not vanish identically on I d . Denote by 

E?(fc) the span of the B-splines M^ s , s € J d {k). If < p < oo, for all kd7L% and all 5 € E«!(fc) 
such that 

5 = a * M k,s, (2-14) 

there is the quasi-norm equivalence 

\\g\\ p x 2-l fc li/l{a s }|| p>fc , (2.15) 



where 



/ \ i/p 

[a s }|| p ,fc := I ^ \a s \ p J 
sGJ r d (fc) 



with the corresponding change when p = 00. 

Let the mixed operator q^, k € be defined in the manner of the definition (|2.6|) by 



% := Il^i ( 2 - 16 ) 



i=l 

We have 

Here and in what follows, for k,k' 6 Z d the inequality k' < k means k\ < hi, i € N[d]. From (I2.17P 
and (I2.10D it is easy to see that a continuous function / has the decomposition 

/ = E ?*(/) ( 2 - 18 ) 



with the convergence in the norm of C(I ). 

From the definition (|2.16p and the refinement equation for the B-spline M, we can represent 
the component functions qk(f) as 

Qk(f)= E c i,sU) M i,s, (2-19) 

where cl are certain coefficient functionals of /, which are defined as follows. We first consider 
the univariate case. We have 

«*(/) = E °M(/) M M- E a k-lAf) M k-l,s- (2.20) 
seJ(fc) sGJ(fc-i) 

If the order r of the B-spline M is even, by using the refinement equation 



M(x) = 2~ r+1 E ( r ^)M(2x-j + r/2), (2.21) 
3=0 
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from (I2.20p we obtain 



q k (f) = Y, c l,s(f) M k,s, 

sg J r (k) 



where 



: = a k,s(f) ~ 4, S (/)> k > °> 
4 iS (/) := 2- r+1 Y *: > 0, a' 0)S (/) := 0. 



(2.22) 
(2.23) 



(m,j)eC(k,s) 



and 



C r (jfe,s) := {(m,j) : 2m + j - r/2 = s, m G J(k - 1), < j < r}, k > 0, C r (0,s) := {0}. 



If the order r of the B-spline M is odd, by using (HOT]) from (123U1) we get (12321) with 

'0, k = Q, 

c fc,s(/) : = < a k,s/2(f), k>0, s even, 

,2- r+1 E( m ,,) e a( M ) GH-l,m(/), fc > 0, a odd, 

where 

C r (A;, s) := {{m,j) : 4m + 2j - r = s, m e J(k - 1), < j < r}, k > 0, C r (0,s) := {0}. 



In the multivariate case, the representation (|2,19p holds true with the c£ which are defined in the 
manner of the definition (|2.7p by 



c k,s(f) — c fci,si(( c fc 2 ,S2(-" C L,Sd(^)))- 



(2.24) 



Let us use the notations: 1 := (1, 1) G R d ; x + := ((x%) + , (x d ) + ) for x G M d ; N d (e) := {s G 
: Si >0, fGe, Si = 0, % £ e} for e C iV[d] (in particular, N d (0) = {0} and N d (7V[d]) = N d ). We 
have N rf (u) n N rf (t>) = if u ^ v, and the following decomposition of Zi: 

Zj= |J N d (e). 

eCJV[d] 

Lemma 2.2 Lei < p < oo and r < min(p, 1). T/ien /or any / G C(I d ) and G N rf (e), there 
holds the inequality 



( 



1/r 



q k (f))\\ P < CJ2 E {2 |s - fc|l/ V(/,2- s ) P } T 



uni/i some constant C depending at most on r,fj,,p,d and ||A||, whenever the sum in the right-hand 
side is finite. 



13 



Proof. From the equality 



= n ( T L~i - t d 



it follows that 



q k = e (-i) H nun^-i = e (-i) M r* 



where /c n is defined by /c" = fcj if % 6 «, and Kf = hi — 1 if i ^ u. Hence, 



ii?fc(/)ii P « E n r * u (/)ii 

uC7V[d] 



(2.25) 



Notice that k u € Z+(u) for some v D e, and < h — h\ < k — k u < 1. Moreover, for s € 

s > if only if s > A;" . Hence, by Lemma 12,11 and properties of mixed modulus smoothness we 

have 

/ \ Vt 



\\T k «(f))\\ P « 



£ | 2 l^li/^(/,2- s ) p } T 

[(«),s>fe' u y 

E {2 |s - fe + |l/ V(/,2~ s ) P } T 



ysez° (v),s>fc« 



£ | 2 l«'-fc|i/V(/,2-( s '-( fc - fc +») p } T 
'ez^(^),s'>fc J 

£ { 2 l^'| 1 /X(/,2- s ) p } T . 



1/t 



The last inequality together with f|2.25j) proves the lemma. □ 

Lemma 2.3 Let < p < oo, < r < min(p, 1), 5 = min(r, r — 1 + 1/p). Then for any f 6 C(I d ) 
and £ Zl(e), i/iere ZioWs i/ie inequality 



u e r (f,2~ k ) p < C 



£ {2^(^) + | 1 ||g s (/))|| p } T 



V 

witt some constant C depending at most on r,fx,p,d and ||A||, whenever the sum in the right-hand 
side is finite. 
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Proof. For simplicity we prove the lemma for e = N[d\, i.e., 1%(e) = 7L%. Let / € C(I d ) and 
k G Zl. From (I2J81) and (pETT]) we obtain 



/ 



1/t 



|AU/)II P < C 



(2.26) 



E ii A u<a/))ii; 

Further, by f)2 . 19|) we get 

AJ-(»(/))= E <j(/)AJiWj). 
jeJ d (s) 

Notice that for any the number of non-zero B-spines in (|2. 19[) is an absolute constant depending 
on r, d only. Thus, we have 

\AU<u(f),*W « E Kj(/)PIAJiWj^)l p , 
jeJ d (s) 

From properties of the B-spline M it is easy to prove the following estimate 



(2.27) 



i d (h) 



\AUM r S)j ,x)\Pdx « 2 



-|s|i-$p|(-log|A|-s)+| 



where we used the abbreviation log \h\ := (log \h\\, ...,log \hd\). Hence, by (|2.27p we obtain 

\A r h (q s (f))\\ p « 2 - 5 K-iogW- s ) + | 12 -N 1 j £ I^C/)!' 1 

\jeJ d (s) 

« 2- 5 K- 1 °sl /l l- s )+l 1 |ks(/)|| P . 



By (gJSD we have 



|AJ-(/)|| P « 



E {2- 5|{ - log|h| - s)+li ii^(/))ii P } r 



From the last inequality we prove the lemma. □ 

For functions / on I d , we introduce the following quasi- norms: 

B 2 (f) := llte(/)}||^(L p) ; 

/ oo ,1/0 

B 3 (f) := E( 2(Q ^ /P)fc |l{ c M(/)}IU) 9 • 
^fc=0 ' 

We will need the following discrete Hardy inequality. Let {a^kez^ and {bk}f.^z d be two 
positive sequences and let for some M > 0, r > 

/>/, < 1/ E (2 <5|(fc_s)+|l a s ) T j • (2.28) 



V s 
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Then for any < f3 < 5, 6 > 0, 

IKMII^ < CM\\{a k }\\ b , (2.29) 

with C = C(f3,6,d). For a proof of this inequality for the univariate case see, e.g, [BJ. In the 
general case it can be proven by induction based on the univariate case. 

Theorem 2.1 Let < p, 9 < oo and 1/p < a < r. Then the hold the following assertions. 

(i) A function f € MB^ e can be represented by the mixed B-spline series 

f = E »(/) = E E c Uf) M ls, (2.30) 

satisfying the convergence condition 

B 2 (f) X S 3 (/) « £(/), 

where the coefficient functionals c r k (/) are explicitly constructed by formula (|2.23p - (|2.24p as 
linear combinations of an absolute constant number of values of f which does not depend on 
neither k, s nor f . 

(ii) If in addition, a < min(r, r — 1 + 1/p), then a continuous function f on I d belongs to the 
Besov space MB^ S if and only if f can be represented by the series (|2.30p . Moreover, the 
Besov quasi-norm B(f) is equivalent to one of the quasi-norms B 2 (f) and B%{f). 

Proof. Since by (I2.15P the quasi-norms B 2 {f) and B%{f) are equivalent, it is enough to prove (i) 
and (ii) for B^(f). Fix a number < r < min(p, 1). 

Assertion (i): For k G Z+, put 

b k :=2^\\q k (f)\\ p , a k := (E {^^(f, 2" fc ) p } T ] 
if k G N d (e). By Lemma [231 we have for 

< c-[e«i x 

s>k 

Then applying the mixed discrete Hardy inequality (I2.28p - (|2.29p with j3 = a — 1/p, gives 
B 3 (f) = \\{b k }\\ b , < C\\{a k }\\ b , x B,(f) X B(f). 

Assertion (ii): Let in addition, a < min(r, r — 1 + 1/p). For k € Zl, put 

1/r 

«fc := lkfc(/)ll P 



6 fc := (E{^(/>2- fc ) P } T ) 
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if k G N d (e). By Lemma 12.31 we have for k G 



l/r 



h < c 



h 5 ^ 



s& 



where 5 = min(r, r — 1 + 1/p)- Then applying the mixed discrete Hardy inequality (|2.28p ~ (|2.29p 
with (3 = a, gives 

B(f) x B x {f) x \\{b k }\\ b , < C\\{a k }U = B 3 (f). 
The assertion (ii) is proven. □ 

Remark From ()2.23j) (|2.24|) we can see that if r is even, for each pair k, s the coefficient c r k s (f) is 
a linear combination of the values f(2~ k (s — j)), and f(2~ k+1 (s' — j)), j G P d (fi), s' G C r (k,s). 
The number of these values does not exceed the fixed number (2/j, + l) d ((r + l) d + 1). If r is odd, 
we can say similarly about the coefficient c r k s (f)- 



3 Sampling recovery 



Recall that the linear operator R m ,m G Z+, is defined for functions on I d in (jl.5j) as follows. 

Rm{f) = £ q k (f) = E E 4, s (/)M fc %. (3.1) 

Lemma 3.1 For functions f on I d , R m defines a linear sampling algorithm of the form (11. ip on 
i/ie grid G d (m). More precisely, 

RmU) = L n {f) = £ /(2- fc J> fcj , 

(fc,s)eG ti (m) 

w/iere 

n := |G d (m)| = E x 2 m m fl! - 1 ; (3.2) 

V'fc j' are explicitly constructed as linear combinations of at most (4// + r + 5) d B-splines M k G 
M r d (m) /or even r, and (12^ + 2r + 13) d B-splines M r k s G M r d (m) /or odd r; M r d (m) := {M£, s , : 
fc' G A(m),s' G J r d (fc')}- 

Proof. Let us prove the lemma for even r. For odd r it can be proven in a similar way. For 
univariate functions the coefficient functionals ak lS (f) can be rewritten as 

ak,s(f) = E Ks-3)fk(2- k j) = E A M (i)/(2- fe j), 
|s-i|<A» ieP(fc,s) 
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where A M (j) := X(s - j) and P(k, s) = P s (jj) := {j G {0, 2 k } : s - j G P(p)} for (i < s < 2 k - 
^k,s(j) is a linear combination of no more than max(r, 2fi + 1) < 2\i + 2 coefficients X(j),j G -P(At), 
for s < /i or s > 2 fc — fi and 



P a (/i)U{0,r-l}, s</i, 
P s ( / u)U{2 fe -r + l,2 fc }, s>2 k -fi. 



P(k,s) C 

Further, for univariate functions we have 

<£,.(/)= E A fe>s (j)/(2-*j) - 2- r+1 £ Q £ A fe _ 1;m (i)/(2- fc (2i)) 

jeP(k,s) (m,u)eC r (k,s) ^ ' jeP(fc-l,m) 

= E A fc , s (j)/(2- fe j), 

j6G(fc,a) 

where G(k,s) := P(k,s) U {2j : j G P(& - l,m), (m,i/) G C(/c,s)}. If j G P(k,s), we have 
|j — s| < max(r, 2/x + 1) < 2/j, + 2. If j G P(k — l,m), (m,u) G C(k,s), we have |2j — s| = 
|2j - 2m - z/ + r/2| < 2| j -m\ + \u- r/2| < 2 max(r, 2/z + 1) + r + 1 < 4/x + r + 5 =: p,. Therefore, 
G(k, s) C Ps(p-), and we can rewrite the coefficient functionals c r k s (f ) in the form 

<£,.(/) = E ^*,-0")/(2- fc j) 

with zero coefficients Xk, s (j) f° r i ^ G(k,s). Therefore, we have 



r 

s 



«&(/) = E <*,.(/ w,« = E E h,sWV- k m r k 

sGj r (fe) s£j r (k) j-s£P(p.) 

= E /( 2 ^') E TfcjW^,. 

je/(fc) s-jeP(p.) 
for certain coefficients 7fe,j (s)- Thus, the univariate %(/) is of the form 

?*(/) = E /( 2 ~ fe i)^' 

where 

V'fcJ := E 7fcj(a)M fcji , 
s-jeP(p.) 

are a linear combination of no more than the absolute number 4^ + r + 5 of B-splines M£ s , and 
the size \I(k)\ is 2 fc . Hence, the multivariate qk(f) is of the form 

?*(/) = E /(2" fc J>fcJ. 
je/ d (fe) 

where 

1=1 
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are a linear combination of no more than the absolute number (4/x + r + 5) d of B-splines MJ, £ 

M!~{m), and the size is 2l fc l 1 . From (|£TJ) we can see that R m (f) is of the form (jl.ip with n 

as in (K2b. □ 



Theorem 3.1 Let < p,q, 9 < oo and 1/p < a < r. Then we have the following upper bound of 
E(m). 



(i) For p > q, 



E{m) <C < 



(ii) For p < q, 



E(m) < 



2~ am , 0<min(g,l), 

2 -Ofn m (d-l)(V9-l/«) i g > m i a (g j 1) ; g < 1 ) 

^ 2 -«m m (d-l) (1-1/0) ^ > m i n (g ) l) ) g > 1. 

' 2 -(a-l/p+l/?)m m (d-l)(l/ g -l/0) +; g < oo, 
2 -(a-l/p)m m (<i-l)(l-l/e) +) (7 = 00. 



Proo/. 

Cose (i): p > q. For an arbitrary / € Ppg, by the representation ()2.30p and (|2. we have 

\\f-Rm(fW q « E 
|fc|i>m 

with any r < min(g, 1). Therefore, if < min(g, 1), then by the inequality \\qk(f) \\g < ||<?fc(/)||p we 
get 

4 ye 

\\f-R m (f)\\ q « ( £ 
\|fc|i>m 

< 2—' f E {2 Q|fc|l |l%(/)ll P } e 

\Jfc|i>m 

< 2- am B 2 (f) < 2~ am 

Further, if > min(o, 1), then 

ii/--R m (/)iif « E iift(/)iif « E {2 Q|fe|i n%(/)ii g } 9 *{2- Q|fe|i r* ! 

fc|i>m |fe|i>m 

where q* = min(o, 1). Putting v = 6/q*, by Holder's inequality and the inequality < 
lkfc(/)llp we obtain 



n/-p m (/)iif « e { 2alk Hqk(f)u g * u \ E {2-° |fc|i r v ' 

\|fc|i>m J Y|fc|i>m y 

< { J B 2 (/)} g *{2~ mn m ( ^ 1)(1/ ' 5, *" 1/e) } 9 * < { 2 -«"i m (d-i)(i/ g *-i/e)| g ^ 



(3.3) 
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This proves the Case (i). 

Case (ii): p < q. We first assume q < oo. For an arbitrary / G B^ e , by the representation (|2.30p 
and Lemma 15.31 we have 

ll/-*m(/)ll| « E {2 {1/P ~ 1/qMl \Mf)\\p} q - 
\k\i>m 

Therefore, if 8 < q, then 

\\f-Rm(f)\\ q « ( E {2 (1/p - 1/9) " £|l ||%(/)ll P } e ) 
\|A|x>m / 
^ 2-{a-l/p+l/q)mg/j:\ ^ 2-( a - 1 /P+ 1 /<2)«\ 

Further, if 6 > </, then 

||/-1W/)||2 « E {2 (1/p ~ 1/9)|fe|l lk(/)U 9 

|fe|i>m 

« ^ {2 Q i fc i i n%(/)ii P r{2- (a - 1/p+1/9)|fc|i r- 

|fe|i>m 

Hence, similarly to (|3.3p . we get 

E q (m) <C {2 _ ( a_1 / p+1 /^ m m^ _1 ^ 1 / IJ_1 / 6, )}' ? . 

When q = oo, the Case (ii) can be proven analogously by use the inequality 

H/-iU/)lloc « E i lkll/p \Mf)\\ P . 

\k\i>m 

□ 

The following theorem for the case a = 1/p can be proven by use of Lemmas 12.21 and 15.31 and 
the inequality (|2.1ip . 

Theorem 3.2 Let < p, q < oo, < 9 < min(p, 1) and a = 1/p < r. Then we have the 
following upper bound of E(m). 



(i) For p > q, 



E{m) < 



' 2 -m/p m (d-l) 



P> 1, 



2 -m/P m (<*-l)/P p < 1. 



(ii) For p < q, 



E(m) < 2" 



m /l m {d-l)/q 
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The following two theorems are a direct corollary of Lemma 13.11 and Theorems 13.11 and 13.21 
Theorem 3.3 Let < p, q, 9 < oo and 1/p < a < r. If m is the largest integer of m such that 

fc€A(m) 

then we have the following upper bound of r n and E{fh). 



(i) For p > q, 

' log d_1 n) a , 6> < min(g, 1), 

(n- 1 log^ 1 n^log^ 1 n) 1 ^- 1 / 9 , 9 > min(g, 1), q < 1, 

_ (n- 1 log^ 1 n) a {\og d - 1 n) 1 - 1 / , 9 > min(g, 1), q > 1. 



r n < £(m) < < 
(ii) For p < q, 



r n « S(m) < 



' (n- 1 log^ 1 n) - 1 ^!/^^- 1 n )(Vff-V»)+ , o < oo, 
(n- 1 log^ 1 n) a - 1 / p (log d " 1 n)( 1_1 /')+ , q = oo. 



Theorem 3.4 Lei < p, < oo, < 9 < min(p, 1) and a = 1/p < r. If m is the largest integer 
of m such that 

then we have the following upper bound of r n and E(fh). 



(i) For p > q, 

r n < E{m) < 

(ii) For p < q, 



n 



- 1 log^ 1 n) l /P log d_1 n, p > 1, 

(n- 1 log^ 1 n) 1 / p (log d - 1 n) 1 ^, p < 1. 

n < £7(m) < (n^log^n^log^n) 1 /". 



From Theorem 13.31 and Lemma l5.ll we obtain the following theorem. 

Theorem 3.5 Let 1 < p,q < oo, < 9 < oo and 1/p < a < r. If in is the largest integer of m 
such that 

i/ien toe /iave i/ie following asymptotic order of r n and E(fh) . 



(i) For p > q and 9 < 1, 

E(m) X r n X log^ 1 n) c 



2 < q < p < oo, 
1 < p = q < oo. 
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(ii) For 1 < p < q < oo, 



E(m) x r n x (n- 1 log d - 1 n) Q - 1 /P+V ?( i og d-i n) (V 9 -i/e) +) J 2 2 < < 1> 

[q < 2. 

4 Interpolant representations and sampling recovery 

We first consider a piecewise constant interpolant representation. Let X[o,i) an d X[o,i] be the 
characteristic functions of the half opened and closed intervals [0,1) and [0,1], respectively. For 
k € Z + and s = 0, 1, 2 k — 1, we define the system of functions N kyS on I by 



N k , s (x) :-- 



X[ ,i)(2 k x-s), 0<s<2 fe -l, 
X[o,i](2 fe x- S ), s = 2 k -l. 

(In particular, iVo.o = X[o,i])- Obviously, we have for k > and s = 0, 1, ...,2 k — 1, 

Nk-i,s = N k)2s + N k)2s +i- 

We let the operator 11^ be defined for functions / on I, for k £ Z+, by 

2 fe -l 

n fc (/) := ^/(2~ fc S )iV M , and II_i(/) = 0. 

s=0 

Clearly, the linear operator life is bounded in L^I), reproduces constant functions and for any 
continuous function /, 

HZ-n-fcC/oiu < ^(/^-^oo, 

and consequently, ||/ — H k {f)\\oo — > 0, when k — > oo. Moreover, for any x G I, IIfc(/, x) = f(2~ k s) 
if x is in either the interval [2 _fc s, 2 _fc (s + 1)) for s = 0, 1..., 2~ fc — 2 or the interval [2~ k s, 1] for 
s = 2~ k — 1, i.e., Ilfc possesses a local property. In particular, interpolates / at the points 

2- k s, s = {0, 1, 2 k - 1}, that is, 

U k (f,2- k s) = f(2- k s), S = 0,l,...,2 fe -1. (4.1) 

Further, we define for k € Z + , 

7r fc (/) := n fc (/)-n fc _i(/). 

From the definition it is easy to check that 

seZi(fc) 

where Zi(0) := {0}, Zi(ife) := {0, 1, 2 k ~ x - 1} for k > 0, 

wLfa) : = ATfc^s+iOc), A: > 0, and ^(x) := iV 0i0 (x), 
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and 

AL(/) := Ai. fc (/,2- fe+1 S ), k>0, andAj j0 (/) := /(O). 
We now can see that every / £ C(I) is represented by the series 

f = E E x lAf)d,s, (4.2) 

fcez + seZi(fc) 

converging in the norm of (I) . 

Next, let us revisit the univariate piecewise linear (nodal) quasi-interpolant for functions on R 
defined in (|2.4[) with M(x) = (1 — |x|)+ (r = 2). Consider the generated from it by the formula 
(|2.5p quasi-interpolant for functions on I 

Q k (f,x) = f(X k s)M k>s (x), (4.3) 

seJ(k) 

and the related quasi-interpolant representation 

/ = E = E E c m(/)^m, ( 4 - 4 ) 

fcez+ fegz + seJ(fc) 

where we recall that J(k) := {seZ: < s < 2 k } is the set of s for which M/% )S do not vanish 
identically on I. From the equality M k>s (2~ k s') = 5 s ^ s i one can see that Qk(f) interpolates s at 
the dyadic points 2~ k s, s € J{k), i.e. 

Q k (f,2~ k s) = f(2- k s), s€J(k). (4.5) 

Because of the interpolation property (|4.ip and (|4.5p . the operators H k and Q k are interpolants. 
Therefore, the representations (14. 2p and (I4.4p are interpolant representations. We will see that the 
interpolant representation (|4.4p coincides with the classical Faber-Schauder series. The univariate 
Faber-Schauder system of functions is defined by 

T:= {<pl t , : s G Z 2 (k), keZ + }, 

where Z 2 (0) : = {0, 1} and Z 2 (k) : = {0, 1, 2 fe ~ 1 - 1} for k > 0, 

Vo.oW : = M o,o(z), Vo,i( x ) : = M o,i( x )> x G 
(an alternative choice is ipo t i(x) := 1), and for k > and s € -Z'(fc) 

^l,s( x ) : = M k,2s+l{x), X G I. 

It is known that J 7 is a basis in C(I). (See |15j for details about the Faber-Schauder system.) 

By a direct computation we have for the component functions q k (f) in the piecewise linear 
quasi-interpolant representation (|4 .41) : 

?*(/) = E x lM)vlM- ( 4 - 6 ) 

sez 2 (fc) 
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where 

Xl s (f) := ~Al k f(2- k+1 S ),k>0, and := /( s ). 

Hence, the interpolant representation (14. 4p can be rewritten as the Faber-Schauder series: 

/ = £ = E E A L(/)^, S ' 

and for any continuous function / on I, 

||/-Q*(/)||oo < W2(/,2-*)oo. 

Put Z r d (fc) := Jill r = 1,2. For fceZf,se Z r d (fc), define 



i=l 



and \ r k s (f) in the manner of the definition (|2.7p by 

^k,s(f) := K 1 ,s 1 ((K2,s 2 (--- :)K k d ,s d (f)))- 

Theorem 4.1 Let r = 1,2, < p, < oo and 1/p < a < r. Then there hold the following 
assertions. 

(i) A function f € MB^ g can be represented by the series 

f = E E a m(/m,*> ( 4 - 7 ) 



converging in the quasi-norm ofMB^g. Moreover, we have 



B*(f) := 









2 (a-l/p)|fe|i / 







< CB(/). 



('mJ /f m addition, r = 2 and a < min(2, 1 + 1/p), then a continuous function f on I d belongs to 
the Besov space MB® e if and only if f can be represented by the series (|4.7p . Moreover, the 
Besov quasi-norm B(f) is equivalent to the discrete quasi-norm B*{f). 

Proof. If r = 2, from the definition and f)4.6j) we can derive that for functions on I d and k G Z+, 
the component function in the interpolant representation (|2.3Up related to the interpolant 

(|4.3p . can be rewritten as 

?*(/) = E x Uf)dJ x )- ( 4 - 8 ) 

sez*(k) 
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Therefore, Theorem 14.11 is the rewritten Theorem 12.11 This does not hold for the case r = 1. 
However, the last case can be proven in a way completely similar to the proof of Theorem 12.11 by 
using the above mentioned properties of the functions ip\ s and operators Hfc. D 

For m G Z + , we have by (|4.8jl 

R r m (f) = R m {f) = Y, E X kM)fl,s- 

For functions / on R r m defines a linear sampling algorithm of the form (|l.ip on the grid G d (m) 
where Gf(m) := {2- fc s : k G A(m), s G I?(k)}, lf(k) := {s G 1 d + : < Si < 2 fc » - 1, i G iV[d]}, 
:= I d (k). More precisely, 

iC(/) = W) = E E f( 2 ~ k iWk,r 



where 



E i*) i * 2r ' 

fegA(m) 



i=i 



and the univariate functions are defined by 



jfe = 0, s = 0, 

fe > 0, s = 2j + 1, 

k > 0, s = 2j, 



and 



fc = 0, 

fe > 0, s = 0, 

k> 0, s = 2j + 1, 

k>0, s = 2 k . 



1 2 

-5<o 

1 f,„2 



From the interpolation properties (|4.1|) and (|4.5|) . the equality ^ s (2 = 5 SjS i one can easily 
verify that R r m (f) interpolates / at the grid G d (m), i.e., 



R r m (f,x)= f{x), xeG d (m). 
Theorem 4.2 Let r = 2, < p,q,0 < oo, and 1/p < a < min(2, 1 + 1/p). Then we have 



(i) For p > q, 



E(m) > 2- am m( d - 1 ^ 1 - 1 /^+ . 
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(ii) For p < q, 



E(m) > 2~( a ~ l l p+1 l q } m rrS d ~ 1 ^ l l q ~ 1 l e ^ + . 



Proof. Put r(m) := {k G N d : |fe|i = m + 1}. Let the half-open d-cube I(k,s) be defined by 
I(k,s) := Y[ i=l d{si2-^- l \{ Si + 1)2~^- 1 )). Notice that J(fc, s) C l d and J(lfe, s) n I(k, s') = for 
s ^ s'. Moreover, if < v < oo, for k G r(m), s £ Z d (k), 



l/u 



I(k,s) 



\VkAx)\ U dx 



-m/v 



with the change to sup when v = oo, and 



E 



<Pk,. 



(4.9) 



(4.10) 



Case (i): For an integer m > 1, we take the functions 



with some fc G r(m), and 



5i .= Ci2 ^ 

sez|(fe) 



fcer(m) sez d (k) 



<Pk,s- 



(4.11) 



(4.12) 



Notice that the right side of (|4.1ip and (|4.12p defines the series (|4.7p of <7j, i = 1, 2. By Theorem 
14.11 and (|4.10p we can choose constants Cj so that gi G -B^g for all m > 1 and z = 1,2. It is easy 
to verify that ^ - R^idi) = 9i i = L 2 - We have by (I4.10P 



S(m) > || 5 i||, > 2- 



if < 1, and 



£(m) > Usall, > || 52 || 9 * > 2- am m( d - 1 )( 1 - 1 / e ) 
if > 1, where q* := min(g, 1). 

Case (ii): Let s(k) G 1% be defined by s(k)i = T.)=i^ j if h > 2, and s(&); = if h = 2 for 
i = 1, d, and T*(m) := {k G r(m) : fcj > 2, i = 1, For an integer m > 2, we take the 

functions 

53 = Ctf-^M^l^ (4.13) 

with some k* G r*(m), and 



04 = c^-y^m-^vo j: < s{k y 

fcer*(m) 



(4.14) 
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Similarly to the functions gi, i = 1, 2, the right side of f|4. 13j) and (|4.14p defines the series (|4.7p of 
gi, i = 3,4, and we can choose constants Cj so that gi £ B® d for all m > 2 and i = 3,4. Obviously, 
gi - R 2 m {gi) = gi, i = 3, 4. We have by <@3§ 

E(m) > \\g 3 \\ q » 2 -^- 1 /P +1 ^ m 

if 9 > q, and 

E(m) > \\g 4 \\ q » 2 - (a - 1 / p+1 M m m (d - 1 ^- 1 ^ 

if 9 < q. □ 

From Theorems 13.11 and 14,21 we obtain 
Theorem 4.3 Zei r = 2, < p,q,9 < oo, and 1/p < a < min(2, 1 + Then we have 



(i) For p > q, 

jV am , 0<mm(g,l), 



E(m) 

(ii) For p < q < oo, 



1 2-am m (d-i)(l-i/e) ) 0>l,g>l. 
£(m) x 2 ~ (a ~ 1/p+1/g)m m (<i ~ 1)(1/<? ~ 1/e)+ . 



Notice that Theorem I4.3f i) has been proven in |22j for the 1 < p = q = 9 < oo. 



5 Appendix 



Lemma 5.1 Let 1 < p, q < oo, < 9 < oo and a > (1/p — l/q)+. Then we have the following 
asymptotic order of X n (B" e ) q . 

(i) For p > q, 



e)q 



[n 




r*" 


"V 


a , 






< 2 < q 


< p < oo 


(n~ 


l \oi 


Jd- 


"V 


a , 






9 < 1, p = 


= q = oo, 


< (n~ 


-Moj 


Jd- 


-y 


a , 






1 < p = g 


< 2, 9 < 


(n~ 


-Moj 


; d- 


~V 




r d -^ n )l/q- 


-1/9 

J 


1 < P = Q 


< 2, 9 > 


(n~ 


Moj 


jd- 




a (l0{ 


; d-l n) l/2- 


-1/0 

5 


9>2. 





(ii) For 1 < p < q < oo, 



'(n- 1 log d - 1 n) Q - 1 /P +1 /<?, 
(n- 1 log^ 1 rc^-VP+i/^logrf- 1 n )(V9-i/0) +) 



2 < p, 2 < 9 < 
q<2. 
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Proof. This lemma was proven in [M], [20] except the cases 6<2<q<p<oo and 9 < 1, p 
q = oo which can be obtained from the asymptotic order [20] 



1. 1 ,,,.,.'/— i ...„» J 1<^<2< (/ <p<oo, 
1 = 1) p = q = oo, 



and the equalities An(W) 9 = An(coW) 9 and coB^q = -B^ max ^ ^, where coVF denotes the convex 
hull of W. □ 

For p = (pi, ...jPd) G (0, oo) d , we defined the mixed integral quasi-norm || • || p for functions on 
I d as follows 



P2/P1 

\f{ X )\^d Xl ) dx 2 --- 



VdlVd 



i/pd 



dx d 



and put 1/p := (1/pi, \jpd). If p, q G (0, oo) rf and p < q, then there holds Nikol'skii's inequality 
for any / G T,f(k), 

< C2l( 1 / p - 1 / q ) fc l 1 ||/|| p (5.1) 



with constant C depending on p, q, d only. This inequality can be proven by a generalization 
of the Jensen's inequality for mixed norms and the following equivalences of the mixed integral 
quasi-norm || • || p . For all k G and all / G S^(/c) of the form (12. 14[) . 



p,A; j 



where 





( 


( 




as} \\p,k '■= 


£ 








1 s d £j(k d ) 


\ 





[PI 



r~; 


Pd/Pd-i\ 

1 ) 







I/Pd 



Lemma 5.2 Let < p < q < 00, 5 = 1/2 — p/(p + q). If k,s G then for any ipk G ^v(^) and 
ip s G Y<f.(s), there holds the inequality 



id 



{x)\ q/2 dx < CA k A s 2^ k ~ sh , 



with some constant C depending at most on p, q, d, where 

A k := (2(^-^)1^11^1^) 



Proof. Put v := (p + q) /p. Then 5 = 1/2 — \ jv and 2 < v < 00. Let v' be given by + = 1. 
Then 1 < u' < 2. Let u, v G (0, oo) d be defined by u := qv/2 and = v if fcj > Sj and V{ = 1/ if 
ki < Si for i = 1, ...,d. Let u' and v' be given by l/u + 1/u' = 1 and l/v + 1/v' = 1, respectively. 
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Notice that v £ (l,oo) d . Applying Holder's inequality for the mixed norm || • || v to the functions 
l^fcl 9 / 2 and l^sl 9 / 2 , we obtain 



\<p k {x)<p 8 (x)\*/ 2 dx < HKI^IIvlH^Klv' = ll^fu /2 ||^||f • 

Since u > pi and u' > pi, by the inequality (|5. 1 j) we have 

\\<Pk\\u < 2K 1 /f-V")fe|i||^|| p; ||^|| u/ < 2 |(i/p-i/«')-|i|| y , 5 || p . 

From (|5.2p and (|5,3p we prove the lemma. □ 

Lemma 5.3 Let < p < q < oo and g € L q be represented by the series 

9 = E E S r(*0- 



1/9 

'■ / ''- L -""l*l^ fc ||| 



Then there holds the inequality 

\\9\U < C I £ ||2CV 

luzf/i some constant C depending at most on p, d, whenever the right side is finite. 

Proof. It is enough to prove the inequality (|5.4j) for g of the form 

9 = 9k, 9k ^ £{?(*), 

fc<m 

for any m £ Zi. 

Put n := [o] + 1. Then < q/n < 1. By Jensen's inequality we have 



k<m 



< i E w*)i 9/n 

, k<m 



9 ( 

V 



k<m 



e-e ni^(*)i 9/n 



consequently, 



By use of the identity 



i JJi^Cx)!'/^. 



l/2(n-l) 



n«i = n« 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



i=i 
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for non- negative numbers a\, a n , we get 

J-= [ fl\ 9k j(x)\*/ n dx = [ X\\g k ,{x)g kJ {x)\il 2 < n ^dx. 
Jld f=i Jld ¥i 

Hence, applying Holder's inequality to n(n — 1) functions in the right side of the last equality, 
Lemma 15.21 and (|5.6|) gives 

j < u([ d \9Ax)g k Ax)\ q/2 dx] 1/n{n 1] < n(^^2- 5|fci ^ |iXl/n(n_1) 

l/2(n-l) ^ 

n 1 
v M?'i'=l j'=l 



l/n(n-in V2(n-1) 



v2/n 



j'=i 



,i'=l 



n4 n n 2 ""****' 11 

i=i Vi=i 



l/n(n-l) 



n n 



1/n 



,i=i i=i 



where A := 6/(n — 1) > 0. Therefore, from (|5.5p and Holder's inequality we obtain 

l/n 

n n 

w ^ e ••• e in^n 2 ~ AMi 

fei^m fc n <m \J=1 i=l 



l/n 



£ II E"E •>• II- 

j=l Vfci^m /c"<m j=l 



•Atfljfej -it'll 



3=1 



(5.7) 



We have 



= E 4, E ■ ■ E E ■■■En 2 " Ml 

fcJ<m k 1 <m ki~ 1 <m ki +1 <m k n <mi=l 



n-1 



E A l E 2 ~ 

ki<m \s<m 



■\5\ki-s\i 



< c E 



fcJ <m 



Using this estimate for £>j, we can continue (|5.7p and finish the estimation of as follows. 



l/n 



< ^E A 



j=l k<m 
k<m 



P {1/p ' 1/qMl g k \\ q p . 
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Thus, the proof of the lemma is completed. □ 

Remark A trigonometric polynomial version of Lemma 15.31 was proven in [25] for 1 < p < q < oo. 
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